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The singlet dynamics which plays a major role in the physics of the spin-1/2 Quantum Heisenberg 
Antiferromagnet (QHAF) on the Kagome lattice can be approximately described by projecting 
onto the nearest-neighbor valence bond (NNVB) singlet subspace. We re-visit here the effective 
Quantum Dimer Model which originates from the latter NNVB-projected Heisenberg model via a 
non-perturbative Rokhsar-Kivelson-like scheme. By using Lanczos exact diagonalisation on a 108- 
site cluster supplemented by a careful symmetry analysis, it is shown that a previously-found 36-site 
Valence Bond Crystal (VBC) in fact competes with a new type of 12-site " resonating-columnar" 
VBC. The exceptionally large degeneracy of the GS multiplets (144 on our 108-site cluster) might 
reflect the proximity of the Z2 dimer liquid. Interestingly, these two VBC "emerge" in different 
topological sectors. Implications for the interpretation of numerical results on the QHAF are outlined. 

PACS numbers: 75.10.Kt, 75.10.Jm, 75.40.Mg 



I. INTRODUCTION 

The spin-1/2 quantum Heisenberg antiferromagnet 
(QHAF) on the kagome lattice (see Fig. 1) is the 
paradigm of frustrated quantum magnetism. Herbert- 
smithite is an excellent experimental realization of such 
a kagome QHAF : Copper atoms carrying spin-1/2 de- 
grees of freedom interacting antiferromagnetically are lo- 
cated on the sites of perfect kagome layers. In addition, 
deviations from an exact QHAF model like inter-layer 
couplings, Dzyaloshinski-Moriya magnetic interactions 2 , 
intrinsic impurities 3 , etc... remain weak. The absence of 
any magnetic ordering at the lowest attainable temper- 
atures 4 is the whole-mark of a non-magnetically ordered 
ground state as it may be realized in the ideal QHAF 
theoretical model (see below). However, it is still diffi- 
cult to say whether the absence of a spin gap as inferred 
from the low-temperature behavior of the susceptibility 5 
can also be interpreted as a generic feature of the perfect 
(theoretical) QHAF or whether it is due to an extrinsic 
perturbation, as one of those mentioned above. 

On the theoretical side, the proliferation of many low- 
energy singlets revealed in early Lanczos Exact Diagonal- 
isations (LED) of small clusters 6 strongly suggests that 
a non-magnetically ordered phase could be stable in the 
Kagome QHAF. Spin liquids (SL) preserving all sym- 
metries, both spin-SU(2) and lattice translation, are the 
most fascinating and challenging candidates. Algebraic 
SL 8 or gapped SL 9 ' 10 could be realized. Recently, Den- 
sity Matrix Renormalization Group (DMRG) results 11 on 
large Kagome strips found strong signals of a Z 2 gapped 
SL. However, in Variational Monte Carlo (VMC) stud- 
ies 12 Z2 SL (whose fermionic formulations were recently 
proposed 13 ) were found to have slightly higher energy 
than the algebraic SL. 

Valence Bond Crystals (VBC) are other exotic candi- 
dates which, in contrast to SL, spontaneously break the 
elementary (3-site) unit cell translation symmetry, hence 
realizing a spontaneous (small) modulation characterized 



by a larger unit cell (named here "supercell" ) . VBC with 
12-site 14 , 18-site 15 , 6-site 16 and 36-site 15 ' 17 ' 18 super- 
cells have all been proposed in the literature as possible 
ground state (GS) of the Kagome QHAF. Very schemat- 
ically, a VBC is often drawn as a frozen hardcore cov- 
ering of the lattice by nearest-neighbor spin singlets, so 
called nearest neighbor valence bonds (NNVB). In real- 
ity, strong quantum fluctuations are expected to severely 
reduce the dimer order parameter. As shown in Fig. 1(a), 
the 36-site VBC shows an hexagonal lattice of resonating 
hexagons. It corresponds to a rather large 2\/3 x 2\/3 
super-cell. The 6-site VBC exhibits a 2 x 1 super-cell 
with columnar dimer order as shown in Fig. 1(b). Both 
VBC are of particular interest in this study and we shall 
refer to them as columnar VBC (CVBC) and hexago- 
nal VBC (HVBC). In VMC studies, a (fermionic) HVBC 
can be stabilized by a very small next NN ferromagnetic 
exchange coupling. 19 The observation of "diamond" pat- 
terns characteristic of the diamond VBC (DVBC) dis- 
played in Fig. 2(b) (deriving from the "parent" VBCo 
of Fig. 2(a)) was observed in recent DMRG 11 studies. 
Also, closely related to VBC , the 12-site VBC of Ref. 14 
shown in Fig. 2(c) displays a triangular lattice of res- 
onating "stars" and shall be denominated as star VBC 
(SVBC). 



II. GENERALIZED QUANTUM DIMER MODEL 

All these theoretical studies show that singlet dy- 
namics plays a major role in the low-energy physics 
of the spin-1/2 Kagome QHAF. One possible route is 
therefore to approximate the QHAF by projecting onto 
the NNVB singlet non-orthogonal subspace. 20 Starting 
from the NNVB-projected Heisenberg model, an effective 
Quantum Dimer Model (QDM) can be constructed via a 
non-perturbative Rokhsar-Kivelson (RK) like scheme. 21 
This procedure bears the enormous advantage to provide 
a mapping to an orthogonal hard-core dimer basis. De- 
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FIG. 1: Kagome lattices. The 108-site clusters delimited by 
the (red) dashed lines are made of 6 x 6 unit cells of 3-site 
up-triangles. Periodic boundaries are used. The four topo- 
logical sectors are determined by the parities of the number 
of dimers (black thick bonds) cut by the three a-a b-b and 
c-c closed loops along the three lattice directions shown (see 
text). For clarity, the hard-core dimers are only shown along 
the three cuts. "Defect" triangles with no dimers are left 
empty, (a) 2\/3 x 2\/3 (36-site) HVBC showing resonating 
hexagons (shaded in green) and stars (shaded in blue) in the 
[1,1,1] topological sector, (b) 2 x 1 (6-site) CVBC with two 
types of up and two types of down triangles (all colored differ- 
ently) in the [0,0,1] topological sector. The red stars indicate 
the inversion (^-rotation) centers. 



tails can be found in Refs. 22,23. Due to the orthogonal 
nature of the dimer basis, the Kagome effective QDM can 
be diagonalized (by the Lanczos algorithm) on clusters up 
to 108 sites (6x6 3-site unit cells) providing clear signa- 
tures of the HVBC and the proximity (in some parameter 
space) of a quantum transition to a Zi dimer liquid 24 (the 
hard-core dimer equivalent of the SL). Note that, strictly 
speaking, in the QDM framework the spinon (i.e. the 
elementary S=l/2 magnetic excitation) gap is pushed to 
infinity although some attempts have been carried out to 
bring spinons back into the picture. 25 

A remarkable and useful feature of the effective QDM 
is that it provides simple and direct access to topological 
properties of the ground state, which is more involved 
in SU(2)-models. On a periodic torus the parity of the 




FIG. 2: VBC with 2x2 supercells in the [1,1,1] topological 
sector of the periodic N = 108 site cluster, (a) VBCo with 
4 types of up (and down) triangles displayed with different 
colors except the "defect" triangles left empty (degeneracy 
g — 8). (b) The resonances of DVBC (colored diamonds) 
restore reflection symmetry w.r.t. the vertical and horizontal 
directions but suppress the 7r/3-rotation symmetry (g — 12). 
(c) SVBC has higher symmetry than VBCo, resonating stars 
restoring reflection symmetry w.r.t. the crystallographic axes 
(5 = 4). 



numbers of dimers cut by two closed lines going along 
two orthogonal directions across the torus is conserved 
by the QDM Hamiltonian, hence defining four indepen- 
dent topological sectors. As clearly seen in a "toy" RK 
QDM defined in Ref. 26, the Z 2 dimer liquid is character- 
ized by topological symmetry breaking with 4 degenerate 
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GS, one in each topological sector. The Z 2 SL also car- 
ries (gapped) excitations which are vortices of the dimer 
field (visons). In a cluster with an even (odd) number 
of sites, visons come in pairs 24 (in odd numbers 27 ). In- 
terestingly, a class of dimer liquid- VBC quantum phase 
transition can be described by a specific field theory 28 
and has been observed numerically by interpolating the 
Hamiltonian between the toy RK QDM and the effective 



Kagome QDM. 24 

In the present work, we shall use topological symme- 
try extensively to re-visit the effective QDM. From this 
refined analysis, we find that the previously-identified 36- 
sitc HVBC competes with new VBC. 

We start here with the effective QDM derived in Ref. 23 
(see Eq. (42) of this paper), 




The open diagrams correspond to kinetic processes 
around closed loops of length 6 (hexagons), 8, 10 or 12 
(stars). In such a diagram, an alternating covering of 
(hard-core) dimers on the loop (i.e. a dimer on every 
second bond) is moved collectively by one step along 
the loop (i.e. it is "flipped"). Each diagram includes 
of course the reverse flipping action. Also, by definition, 
each diagram is meant foanalysr all (non-equal) loops of 
the lattices with the same topology i.e. obtained under all 
possible translations, reflections and rotations (i.e. under 
all elements of the space group of the lattice) . The sec- 
ond kind of diagrams (colored ones) are diagonal terms 
which simply "count" the number of "flippable" loops 
of the first kind. More details can be found e.g. in 
Refs. 22,23. Note that, in contrast to Ref. 24, Eq. (1) 
includes infinite order resummation of all the processes 
(loops) restricted to a single hexagon. This however only 
leads to very small changes of the values of the coeffi- 
cients of H e g compared to Ref. 24. Note also that the 
stars have vanishing amplitudes at all orders. 



III. TOPOLOGICAL PROPERTIES AND 
SYMMETRY ANALYSIS 

The generalized QDM is investigated by Lanczos ex- 
act diagonalization of the finite 108-site cluster of the 
Kagome lattice made of 6 x 6 unit cells of 3-site up- 
triangles and delimited by (red) dashed lines in Figs. 1,2. 
Periodic boundary conditions are used in the horizon- 
tal a-a and b-b (at 60°-degrees) directions so that the 
system is topologically equivalent to a torus, providing 
global invariance under the translation group T . On such 
a geometry, the generalized QDM hamiltonian bears in- 
teresting topological properties: its Hilbert space of (or- 
thogonal) NNVB hardcore coverings can be split in four 
disconnected topological sectors which are conserved by 
the QDM hamiltonian. To determine the topological sec- 
tor of a given dimer configuration, one simply needs to 



count the number of dimers (black thick bonds) cut by 
the three a-a b-b and c-c closed loops along the three 
lattice directions shown in Figs. 1,2. We shall label these 
sectors as [a,/?,7] where a, f3 and 7 are or 1 depending 
whether the numbers of dimers cut by the correspond- 
ing loops are even or odd, respectively. Every cluster of 
P x P unit cells (with N — 3P 2 sites) contains a "sym- 
metric" topological sector [a, a, a] with a — or a = 1 
for P = 4p or P = 4p + 2, respectively. Flipping the 
dimers along a close loop which goes around the torus, 
e.g. along the a-a direction, switches the parity of the 
number of dimers cut in the two other directions, e.g. 
along b-b and c-c. For the 108-site cluster of interest here, 
one can then build the 4 topological sectors T = [1, 1, 1], 
Ti = [1,0,0], T 2 = [0,1,0] and T 3 = [0,0,1]. All con- 
figurations of each topological sector can be obtained 
separately by constructing an initial dimer configuration 
belonging to the targeted sector and then by applying 
recursively loop- flipping operators of the type of Eq. (1) 
to generate (numerically) new configurations until the 
Hilbert space is complete. 

In order to block-diagonalize the Hamiltonian matrix 
in each topological sector, one can make use of all space 
group symmetries leaving each sector globally invariant. 
The corresponding point group Q{T a ) depends on the 
topological sector T a . Note that we use here the center 
of the hexagon as the common origin for all point group 
elements. The full point group Q — C§ v of the infinite 
lattice can be used for the "fully symmetric" To = [1, 1, 1] 
topological sector. In contrast, the 2-7r/3-rotation is lost 
in the T x = [1, 0, 0], T 2 = [0, 1, 0], and T 3 = [0, 0, 1] topo- 
logical sectors so that the invariant point group is reduced 
to Q = C2V ■ Each of these 3 topological sectors and their 
two corresponding (orthogonal) reflection symmetry axes 
are related by 27r/3- and 4-7r/3-rotations. 

In practice, one can specify a fixed momentum K in 
the Brillouin zones (BZ) of Fig. 3 and diagonalize the 
Hamiltonian separately in all irreducible representations 
(IRREP) of its little group Gk (which is defined by all the 
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TABLE I: Little groups Gyl for all K momenta of Fig. 3(a) 
and Fig. 3(b) corresponding to the To = [1,1,1] and Tj = 
[0, 0, 1] topological sectors, respectively. 



rections (consistently with Ref. 7) while, for the center 
r*, it refers to the reflection along the axis perpendicu- 
lar to the direction c (i.e. along the momentum A± in 
reciprocal space). 



IV. SPONTANEOUS SYMMETRY BREAKING 
AND DEGENERACIES OF VBC 



FIG. 3: 36 available momenta of the 108-site cluster (6x6 
up-triangles) in the first Brillouin zone. Topological sectors 
have to be distinguished (as done on graph) because of their 
different point groups (see text). Equivalent momenta are 
represented by the same (filled) symbols. Translation by a 
reciprocal lattice vector (thin lines) of the momenta at the 
border of the 1st BZ are shown by open symbols (correspond 
to the same momenta). The dashed lines in (b) represent the 
two (orthogonal) reflection symmetry axes of the C 2v group. 



elements of Q(T a ) leaving K invariant in the BZ). All lit- 
tle groups are listed in Table I for the relevant K-points 
of the TV = 108 cluster. The size of the Hilbert space 
(i.e. the # of orthogonal states) can then be reduced 
from the full size 2 N / 3+1 ~ 137 x 10 9 by (i) a factor 4 
for each topological sector, (ii) a factor N/3 = 36 by us- 
ing translation symmetry and (iii) a factor card{^x} (see 
Table I) by using point symmetry. Hence, after proper 
construction of the new "symmetrized" basis states in 
each IRREP, the linear sizes of the corresponding Hamil- 
tonian "blocks" range from ~ 80 x 10 6 (<?k = C§ v ) to 
~ 960 x 10 6 (£ K = Id). In the following, each IRREP will 
be labelled by its momentum/topological sector accord- 
ing to Fig. 3 (different combination of letters, subscripts 
and superscripts are used to distinguish all of them) and 
by the characters (ra, r2, a) which refer to 2ir/3 rotations, 
parity under inversion and reflection about the momen- 
tum direction, respectively, whenever such symmetries 
are relevant (depends on Gk)- For the center of the BZ 
r, a refers to the reflection along one of the lattice di- 



The spectrum of a system which spontaneously breaks 
some lattice symmetries at zero temperature shows some 
ground-state degeneracies. Moreover, these degeneracies 
occur between levels belonging to different IRREP of the 
lattice symmetry group. The simplest example is that 
of a Z2 symmetry breaking, where a state which in the 
"even" (topological) sector is degenerate with a state in 
the "odd" sector. In a finite size-system these degen- 
eracies are only approximate, but they are expected to 
become exponentially precise when increasing the sys- 
tem size. In turn, the presence of degeneracies is al- 
most always the indication of some spontaneous symme- 
try breaking, and this type of signature has been widely 
used in exact diagonalization studies to detect ordered 
phases. 7 

If one knows the VBC pattern/symmetry, it is rela- 
tively easy to compute the IRREP which should appear 
in the ground-state multiplet. Let us denote by |1), • • • \d) 
the states obtained by putting the VBC pattern on the 
lattice in all the possible positions/orientations. For a 
sufficiently large lattice they can be considered as orthog- 
onal to each other, and they transform into each other un- 
der the action of the lattice symmetry group G — Q 
for any symmetry g € G and for any state there exist 
another state j = f(g,i) such that g\i) = \j). We as- 
sume in particular that the later equation holds without 
any phase factor. Such a ground-state multiplet defines 
nothing but a (reducible) representation of dimension d 
of the group G. Determining the associated ground-state 
quantum numbers amounts to decompose this represen- 
tation over the IRREPs of G, which is a standard task in 
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group theory. If we denote by p some IRREP of G, the 
number n p of ground-states which should appear in the 
symmetry sector p is given by a character formula: 
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W\ 



^2x P (g 1 )Tr|i>,--|d>(ff) 
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where Xp(5 X ) i s the character of p on the group ele- 
ment <7 _1 , \G\ is the total number of symmetries, and 
Trix),... \d){g) is the number of VBC states which are in- 
variant under g: g\i) = The later quantity only de- 
pends on the symmetries of the VBC and not on the full 
wave-function. It can therefore be computed using "sim- 
plified" states with the correct symmetries. It can some- 
times be reduced to a single (non-resonant) covering, or a 
linear combination of 2" ceU coverings for resonating states 
with n co n super cells (n cc n = N/N). In any case, from 
the group characters Xp an< i the symmetries of the VBC, 
the multiplicities n p can be systematically computed and 
compared to the spectrum obtained in LED. 
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TABLE II: Quantum numbers of the GS multiplet for VBC 
with m x n super-cells of N = 3mn sites. The letters refer to 
the momenta of Fig. 3 (with their multiplicities in brackets) in 
the two topological sectors [1,1,1] and [0,0,1]. The character 
related to 27r/3 rotations (r 3 = l,j,j 2 ), parity under inversion 
(r 2 = +, — ) and reflection about the momentum direction 
(<t = +, — ) are denoted as (r3,r 2 ,<r). "±" means that both 
even and odd GS are present and "x" stands for "symmetry 
not relevant" (e.g. 27r/3 rotation and hence r?, is irrelevant 
in the [0,0,1] topological sector). VBC , DVBC and SVBC 
are 8-fold, 12-fold and 4-fold degenerate, respectively. CVBC, 
rVBCi, rVBC'i and rVBC3 have the same degeneracy g = 
8 in the [0,0,1] sector i.e. g — 8 x 3 = 24 including all 3 
topological sectors while rVBC 2 has g = 16 x 3 = 48. The 
last line corresponds to the Z 2 dimer liquid (g — 4). 

We have computed these multiplicities for a few VBC 
with a 2 x 2 or 2 x 1 supcrccll on the 108-site cluster. 
These multiplicities were obtained in a fixed topological 
sector, in order to compare to the LED results (next sec- 
tion). To do so, the lattice symmetry group was simply 
restricted to the elements which conserve the topological 



sector, i.e excluding 27r/3 rotations in the [0, 0, 1], [0, 1, 0] 
and [1,0,0] sectors. The results are shown in Table II. 
The VBCs named CVBC is the simplest, since, as far as 
the symmetries are concerned, it can be reduced to a sin- 
gle dimer covering. CVBC is displayed in Fig. 1(b), and 
has first been proposed in Ref. 16. It has a 2 x 1 supercell 
and a degeneracy equal to 8 (for each of the 3 topological 
sectors, so 24 in total). The VBCo of Fig. 2(a) can also 
be reduced to a single dimer covering but, in contrast to 
the CVBC, belongs to the [1,1,1] topological sector. Its 
degeneracy is 24. The other VBC we considered involve 
some resonance between different dimer configurations. 
The rVBC2 (see Fig. 4(a)) can be constructed from the 
CVBC configuration by applying some resonance opera- 
tors a x (notation of Ref. 26) on a quarter of the hexagons 
(colored in blue in Fig. 4(a)), to shape a 2 x 2 supercell: 

|rVBC 2 )= Yl (l + a^)|CVBC). 

h=l,—,9 

Doing so simultaneously enlarges the supercell and dou- 
bles the degeneracy (16). If the coefficient a is set to 
one, one restores a reflection symmetry (cr = +) and 
we get another VBC, named |rVBCi) and with degen- 
eracy 8. By applying the resonance operators on the 
other hexagon-sublattice (colored in orange in Fig. 4(b)), 
one obtains another resonating VBC, IrVBC^) shown 
in Fig. 4(b), with the same degeneracy (8) but differ- 
ent quantum numbers (the symmetry axis is now that of 
the wave vector B* i.e. along the c direction). Another 
resonating crystal, |rVBC3) shown in Fig. 4(c), with no 
reflection symmetry (in contrast to the previous ones) 
but with inversion symmetry w.r.t the centers of all the 
hexagons can also be constructed in the [0,0,1] topolog- 
ical sector but using longer range dimers. We believe 
there is a similar VBC (that we name VBC*) which bears 
exactly the same symmetry property and which could 
be represented solely in terms of (hardcore) NN dimers, 
probably requiring a more involved resonance structure. 
Finally, the SVBC discussed in Ref. 14 is made of res- 
onating "stars" (twelve-dimer loops) and also has a 2 x 2 
supercell. Similarly, the DVBC is made of resonating 
diamonds (eight-dimer loops). But unlike the previous 
examples, these two crystals belong to the [1,1,1] sector. 

It appears that the trivial IRREP (L(+, +, +) or 
L* ( x , +, +)) appears exactly once in all the VBC we have 
listed so far. This is in fact a general property of such 
kind of broken symmetry phase. Indeed, the multiplicity 
of L is just a special case of Eq. 2 where the character \r 
is constant and equal to unity. We thus have 

n r = t^t Tr |i>,-|d>(3)- 

sec? 



|G| 



However Trin ...i^) (<?) is nothing but the number of states 
which are invariant under the symmetry g. So nr can be 
written using the stabilizers Stab(i) or each state \i): 



n r = 
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V. NUMERICAL RESULTS 






FIG. 4: Resonant VBC with 2x2 supercells in the [0,0,1] 
topological sector of the periodic N = 108 site cluster. Res- 
onating 10-dimer (a,c) and 8-dimer (b) loops (see text) are 
represented by shaded areas, (a) When a = 1 (see text) a 
symmetry axis perpendicular to the c direction is restored 
(light and dark green triangles become equivalent), (b) The 
pattern exhibits a reflection symmetry w.r.t. the c direction, 
(c) Inversion symmetry w.r.t. the centers of all the hexagons 
is present. Dimers beyond NN have been used. 



We have used the Lanczos Exact Diagonalization 
(LED) technique to obtain the low energy spectra of the 
QDM Hamiltonian in each IRREP. Computations have 
been performed on the Altix SGI cluster at CALMIP 
(Toulouse), in particular on the Altix UV node provid- 
ing 96 CPUs allocating 1024 Gb of fast access RAM. 
Typically, only 200 Lanczos iterations were sufficient to 
get at least 5 or 6 of the lowest eigenenergies (in each 
IRREP) with full machine (double) precision. Fig. 5(a) 
shows the obtained low-energy spectra at the most rele- 
vant momenta and in the two topological sectors [1,1,1] 
and [0,0,1]. For each of these sectors, we can see a group 
of quasi-degenerate levels named "ground state mani- 
fold" and a dense accumulation of levels at higher en- 
ergy. This clearly reflects the presence of two continua 
of excited states above two (slightly different) thresholds 
separated from the GS manifolds by two gaps. Note that 
the energy spectra at the other momenta (D, E and F 
in Fig. 3) are not shown on this picture since they only 
contain excited states beyond the gaps and no level be- 
longing to the ground state manifold. Fig. 5(b) shows a 
zoom in of the ground state manifold which reveals two 
subgroups corresponding to each topological sector. 
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FIG. 5: (a) Low energy spectrum of the 108-site cluster in 
the T = [1,1,1] and T 3 = [0,0,1] (degenerate with [0,1,0] 
and [1, 0, 0]) topological sectors. Different symbols distinguish 
different momenta and topological sectors (see Fig. 3(a) & 
(b)) and IRREPs labeled as " (r 3 , r 2 , cr)" (see text). Only a few 
levels corresponding to the on-sets of the two continua (shaded 
areas) are shown, (b) Zoom in of the (quasi-) ground states. 
An Ising degeneracy (see text) is reflected by the (almost) 
perfect equal spacings between the HVBC ground states. 



But since all the stabilizers have the same number of 
elements, |Stab(i)| = \G\/d, we just find nr = 1. In 
Sec. VI this remark will allow us to show that a particular 
multiplet of low-lying states cannot be interpreted as a 
single conventional VBC. 



In agreement with Ref. 24, the quantum numbers of 
the GS levels in the "symmetric" [1,1,1] topological sector 
correspond exactly to those predicted (see Ref. 7 for de- 
tails) for the 36-site HCVB schematically shown in Fig. 1. 
In addition, we find here some very remarkable (Ising- 
like) fine structure with equally spaced sub-subgroups of 
levels which shall be discussed later on. The numerical 
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values of these eigenenergies reported in Table III reveal 
almost exact degeneracies within less than 10~ 5 (in units 
of the exchange constant J). 
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energy levels 


degen.* 


r(i,+,-) 


1 


79 548 096 


2 


-3.65721441 
-3.649 849 41 


2 


r(i,+,+) 


1 


79 548 096 


2 


-3.653 462 67 
-3.646 352 04 


2 


A(x,+,-) 


3 


238 642176 


2 


-3.657 214 51 
-3.649 849 60 


6 


A(x,+,+) 


3 


238 642176 


2 


-3.653 462 48 
-3.646 35192 


6 


B(l, x,-) 


2 


159 073 536 


2 


-3.657 213 37 
-3.649 849 46 


4 


5(1, x,+) 


2 


159 073 536 


2 


-3.653 462 55 
-3.646 352 03 


4 


C(x,x,-) 


6 


477218 560 


2 


-3.657 213 33 
-3.649 849 58 


12 


C(x, x,+) 


6 


477 218 560 


2 


-3.653 462 56 
-3.646 35192 


12 



TABLE III: Quasi-degenerate GS classified according to the 
IRREP of the Cq v <S> T space group of the [1,1,1] topological 
sector of the 108 cluster. ^The Lanczos algorithm does not 
(easily) provide the exact degeneracy of each energy level. 
'An extra (exact) degeneracy of 2 is expected for each level 
(see text). 

The investigation of the [0,0,1] and the two other equiv- 
alent topological sectors was overlooked in Ref. 24. In- 
terestingly, Fig. 5(b) reveals a new subgroup of energy 
levels within - 6 x 10~ 3 J from the [1,1,1] GS sub group. 
Their energies are listed in Table IV showing again al- 
most perfect degeneracies within sub-subgroups of levels. 
The quantum numbers of these states can be straightfor- 
wardly compared to the expected quantum numbers for 
the VBC shown in Table II. Surprisingly, one finds that 
the numerical results do not correspond to any single 
VBC. However, they correspond exactly to the reunion 
of the GS multiplets of CVBC and VBC*. 



VI. DISCUSSIONS 

Our results clearly show a close competition between 
two (or even three) types of VBCs. We have checked that 
very small changes of the Hamiltonian parameters can 
easily reverse the relative position of the two subgroup of 
levels. Therefore it is not clear which of these VBC will 
be stabilized in the thermodynamic limit. 

Let us first discuss the remarkable fine structure of the 
GS manifold in the [1,1,1] topological sector. It reflects 
the hidden Ising variable 18 associated to the two possi- 



IRREP # of # of # of GS total 



K states levels t energy levels degen. 





1 


238 613 056 

iUU VJ-LO \_JiJ\J 


2 


— 3 639 546 99 

— 3 639 460 27 


2 


r*(x,+,-) 


1 


238 613 056 


2 


-3.639 546 42 
-3.639 458 83 


2 




1 


238 613 056 


1 


— 3 639 497 82 


1 




1 


238 613 056 


1 


—3 639 495 81 


1 


Ax (x,+,+) 


1 


238 605 760 


1 


-3.639 509 03 


1 


Ai(x,+,-) 


1 


238 605 760 


1 


-3.639 509 00 


1 


A 2 (x,+,x) 


2 


477 218 816 


3 


-3.639 508 79 
-3.639 508 60 
-3.639 49719 


6 


A 2 (x,-, x) 


2 


477218 816 


1 


-3.639 496 79 


2 



TABLE IV: Quasi-degenerate GS classified according to the 
IRREP of the C 2v ® T space group of the [0, 0, 1] topological 
sector of the 108 cluster. ^The lanczos algorithm does not 
provide (easily) the exact degeneracy of each energy level. 



ble chiralities, let's say f and J,, of each of the resonat- 
ing stars (or pinwheels) of the HVBC. On our 108-site 
cluster 3 such resonating stars can be fitted leading to 
2 3 = 8 possible configurations. The numerical spectrum 
can be explained by assuming a (small) effective "mag- 
netic" Ising field h which naturally provides splitting by 
energies — 3/i, —h, h and 3/i. Although this structure is 
relatively robust and does not require real "fine tuning" 
of the effective model, the fact that J12 = is impor- 
tant. Indeed, when J12 is added by hand, the numerical 
results can be simply interpreted by assuming a (strong) 
linear dependance of the effective Ising field with J12, 
h = ho + AJ12, and h changes sign for a tiny (posi- 
tive) value of J12 ~ 0.001. Similarly, a small value of 
J 12 = ±0.01 increases the splitting within the GS mani- 
fold by a factor ~ 5. The behavior in the thermodynamic 
limit (for Jyi — 0) depends crucially whether the effec- 
tive ho is exactly zero or retains a very small but still 
finite value. In the latter (most probable) case, the ab- 
solute HVBC GS should exhibit a uniform chirality of 
all resonating stars (one of the two ferromagnetic Ising 
configurations) and the VBC gap should be filled by a 
discrete set of energy levels corresponding to the excita- 
tions of a finite number of resonating stars. In contrast, 
in the unlikely case where ho = (due to fine-tuning) the 
HVBC GS would retain an infinite degeneracy. 

Identifying the exact nature of the phase in the 2\ , T 2 
and T3 topological sectors is more subtle. In fact, be- 
cause of two (GS) energy levels in the fully symmetric 
r* ( x , +, +) IRREP this phase can not be described as a 
unique VBC. From a direct comparison of the quantum 
numbers of the states in the GS manifold with those listed 
in Table II, our numerical data can indeed be interpreted 
as the simultaneous occurrence of two very closely com- 
peting VBC, the columnar CVBC and the 2 x 2 VBC*. 
Such a description however raises a number of concerns. 
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First, the exact nature of VBC* is not fully known: it 
is similar to SVBC in terms of symmetry properties but 
belongs to a different topological sector. A naive rep- 
resentation in terms of a simple NN dimer covering is 
therefore impossible and should involve subtle resonances 
like rVBCl or rVBC2. Secondly, we have checked that 
the almost exact degeneracy between CVBC and VBC* 
is robust and survives under small changes of the QDM 
Hamiltonian parameters so that it does not correspond 
to any "fine tunning" . This unconventional feature might 
have different origins: (i) CVBC and VBC* are truly dif- 
ferent but are related by some hidden symmetry of the 
original projected Heisenberg model or of its approximate 
effective model (f); (ii) The relevant phase still corre- 
sponds to a unique VBC of smaller degeneracy (3x8 = 24 
instead of 3 x 16 = 48) but admits some form of zero- 
energy excitations. It is interesting to note that (1) only 
contains loops of different lengths but always involving a 
single hexagon. It might be that smaller terms involving 
longer loops (and hence at least two hexagons) could lift 
this degeneracy, (iii) The observed degeneracy is the sig- 
nature of an ordered phase which is more complex than 
a VBC. We could think, for instance, of some uncon- 
ventional VBC with broken time-reversal symmetry. Al- 
though it is not clear to us how to construct precisely 
such a state, it may have nr = 2. The tight competition 
between several VBC in addition to the existence of a 
Z2 dimer liquid stable very nearby in parameter space 24 
are clear signatures of the highly frustrated nature of the 
generalised QDM model, probably also generic in several 



highly frustrated QHAF. In fact, the observation in re- 
cent DMRG studies of the kagome QHAF 11 of signals of 
the formation of crystalline diamond patterns together 
with finger prints of a Z2 spin liquid suggest that the 
QHAF and the effective QDM might have strong similari- 
ties. In any case, our results reveal an exceptionally large 
quasi-degeneracy of the GS multiplets (96 + 48 = 144 on 
our 108-site cluster !) which might reflect the proxim- 
ity of the Z2 dimer liquid. 24 If the same were to exist 
in a frustrated QHAF, we believe the related LED and 
DMRG results might become quite difficult to analyse. 
Incidentally, our study also shows that the existence of 
quasi-degenerate GS in different topological sectors (on 
the cylinder or on the torus) is not sufficient to com- 
pletely identify a Z 2 dimer (or spin) liquid, our result 
providing a counter-example. 
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